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Heavy hadrons on an anisotropic lattice 
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Results from simulations of quarkonia and heavy-light mesons on an anisotropic lattice are presented. The 
improved quark action and action-parameter tuning used in this study are discussed. 
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1. DESIGNING A 3+1 ACTION 

A fermion action, which is useful for simula- 
tions at extreme anisotropies and has discretisa- 
tion errors of 0(a s a, a%, a\ ), is constructed. The 
details of this action are in Refs |3I4| . Our start- 
ing point is the usual isotropic rotation, used in 
the Sheikholeslami-Wohlert (SW) and D234 ac- 
tions 
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For a fixed Wilson parameter, r (usually equal to 
one) doublers may reappear as a t is made small, 
which hinders simulations with a t <C a s . The so- 
lution proposed here is to treat the temporal and 
spatial directions differently. The temporal ro- 
tations are retained since such one-hop improve- 
ment terms ensure the transfer matrix is posi- 
tive definite and no "ghost" states can propagate. 
Thus Eqs.Hare modified to read 
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and once the covariant temporal derivative is dis- 
cretised the temporal doublers are removed. The 
spatial doublers remain and are removed through 
the addition of an (irrelevant) higher-dimensional 
operator, a Df term, to the Dirac operator. The 
resulting improved action can then be written 
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where M r = [i r jiDi + 7o-Dt + [i r m and /i r 
; I ■ ^ratm). The Wilson-like parameter, s is cho- 
sen such that the doublers receive a sufficiently 
large mass. A lattice discretisation scheme is 
now quite simple. An improved discretisation is 
only needed for the 7j-D, term in Eq. [3] since the 
simplest discretisation would lead to 0(a 2 s ) er- 
rors. Including the gauge fields and the mean- 
link improvement coefficients, u s and Ut the lat- 
tice fermion matrix reads 
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Choosing r = 1 and s = 1/8 yields the sD34 
action proposed in Ref 1.. The authors investi- 
gated the radiative corrections to this action and 
at one-loop in perturbation theory find no con- 
tribution from 0(a s a s m q ) terms, which would 
spoil simulations with heavy quarks. This is a key 
advantage of this action over anisotropic SW- or 
D234-type actions where such terms may appear, 
depending on the choice of the parameter, r 
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2. RESULTS 

Our quenched simulations use an improved 
gluon action designed for precision glueball sim- 
ulations on anisotropic lattices We omit the 
temporal rotations in this exploratory study, leav- 
ing a leading discretisation error of O(at). Fur- 
ther details of the simulation and parameter val- 
ues are in Table ^ A range of masses is investi- 



# gauge configs. 


100 


Volume 


10 3 x 120 


a s 


0.21fm 


a s /r 


0.432 


£ = a s /a t 


6t 



Table 1 

Lattice details. ^The gauge anisotropy was tuned 
nonperturbatively, from the static quark poten- 
tial. 

gated from a t m q = —0.04, 2 close to the strange 
quark mass, to a t m q = 1.5. Data were accummu- 
lated at spatial momenta (0, 0, 0), (1, 0, 0), (1, 1, 0) 
and (1,1,1) in units of 27r/a s L, averaging over 
equivalent momenta. Degenerate and nondegen- 
erate combinations of quark propagators are con- 
sidered. The nondegenerate combination is made 
with the lightest quark, and each of the heavier 
quarks. These mesons are denoted heavy-light 
while the degenerate combination is heavy-heavy. 
Effective masses for both pseudoscalar (PS) and 
vector (V) particles were determined from single 
cosh fits. For all masses and momenta good fits 
to long plateux are observed. Fits to the lightest 
and heaviest PS masses are shown in Fig. ^ 

2.1. Mass dependence of £ 

Quantum fluctuations rcnormalisc the 
anisotropy so the input parameter, £ in the action 
is not the same as the anisotropy determined from 
a physical observable in the simulation. The pa- 
rameter in the action must be "tuned" such that 
the measured anisotropy takes its required value. 
The aim of this study was to investigate the 
mass dependence of the renormaliscd anisotropy, 
examining the extent to which the tuning must 



2 This corresponds to a positive quark mass since Wilson- 
type actions have an additive mass renormalisation. 
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Figure 1. Effective mass plots atm q = —0.04, 1.0 
respectively. For atm q — — 0.04,ti p = (1, 1, 1) the 
signal becomes very noisy at t > 50. 

be repeated as atm q is varied. Both the preci- 
sion of the determination and its deviation from 
the required anisotropy are of interest. We also 
examine the difference between the anisotropy 
determined from PS and V mesons. 

£ was tuned at the lightest mass by measuring 
the slope of the energy-momentum dispersion re- 
lation. A value of six from the dispersion relation 
requires £ = 6.17 in the action. The anisotropy 
was then computed, for a range of masses from 
their dispersion relations. Fig. [2] is a represen- 
tative sample of the dispersion relations for PS 
mesons. For clarity, E 2 (p) — E 2 (0) is plotted as 
a function of (a s p) 2 . In general, very good linear 
(relativistic) behaviour is observed. The mass de- 
pendence of the measured anisotropy is shown in 
Fig. |21 for PS mesons. Similar mass-dependence, 
although with larger statistical errors, is seen for 
the vectors. The plots show that up to atm q ~ 0.3 
the heavy-heavy and heavy-light determinations 
are in agreement. In addition the values are con- 
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Figure 2. Dispersion relations for heavy- heavy 
(left) and heavy-light (right) PS mesons. The 
symbols o, □, O, A are atm q = 0.1, 0.2, 0.5, 1.0 re- 
spectively. 

sistent, within errors, with the target anisotropy 
and are determined to « 3% or better accuracy. 
The charm quark mass on this lattice is close to 
atvn q = 0.2, indicating that charm physics is fea- 
sible computationally and requires little parame- 
ter tuning. 

For atm q > 0.5, the heavy-heavy and heavy- 
light dispersion relations do not give a consistent 
value of the anisotropy. This dependence of the 
measured anisotropy on the input parameter £ in 
the action was studied at fixed quark mass. The 
results, for both PS and V particles are shown 
in Fig 0J The value of £, determined from the 
heavy-heavy dispersion relation, moves closer to 
its target value of six while £ determined from the 
heavy-light physics moves away from this value. 
It is also interesting to note the agreement be- 
tween determinations of £ from pseudoscalar and 
vector particles. The tuning, described above was 
carried out for pions and it is reassuring that al- 
though the vector particles have larger statistical 
errors they nevertheless yield a consistent pattern 
for the mass dependence of £. 

3. FUTURE WORK 

A paper with more details and additional quark 
masses is in preparation We plan to investi- 
gate the feasibilty of precision 6-physics using the 
anisotropic lattice described here. 
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Figure 3. The mass dependence of the anisotropy 
for PS mesons as a function of the heavy-heavy 
meson mass in lattice units. 




Figure 4. The dependence of the renormalised 
anisotropy on the input parameter £ , for PS (left) 
and V mesons (right), o and O are heavy- heavy 
and heavy-light respectively. 
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